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We  believe  significant  progress  has  been  acheived  on  a  number  of  areas 
described  in  our  proposal.  We  have  (1)  successfully  developed  and  applied  a 
theory  to  understand  the  deep  to  shallow  transition  of  the  bandgap  energy 
levels  of  large  chalcogen  complexes  in£S p,  (2)  have  developed  an  ab-initio 
tight-binding-like  electronic  structure  method  for  solids,  (3)  have  made  great 
progress  in  understanding  the  energetics  of  interstitial  impurities  in  compound 
semiconductors,  (4)  have  developed  a  theory  to  predict  equilibrium 
concentrations  of  intrinsic  and  extrinsic  defects  in  semiconductors,  (5)  have 
performed  the  first  molecular  dynamics  simulations  of  a  compound  semiconductor 
surface,  and  (6)  have  studied  the  hydrogen  halide  molecular  crystals  under 


pressure.  A  brief  summary  of  these  topics  will  now  be  given. 


r 


A-  theory  for  the  bandgap  levels  of  large  defect  complexes  has 
been  developed  and  applied  to  large  (up  to  five  atoms)  chalcogen  S,  Se,  or  Te 
complexes  in  Si.  It  is  found  that  the  number  of  bandgap  levels  is  equal  to  the 
number  of  impurities  in  the  cluster.  Most  of  these  levels  are  however 
electrically  inert,  except  for  the  the  upper  pair  of  levels.  We  find  that  the 
wavefunction  character  for  these  "active”  levels  is  very  different  to  that  of 
the  isolated  impurity.  The  theory  predicts  that  the  "active"  uppermost  level 
continuously  changes  from  being  deep-level-like  to  more  effective-mass-like  as 


the  size  of  the  defect  complex  incj 


;.  A  number  of  intra-center  transitions 


between  localized  bound  states  have  been  predicted.  The  theory  explains 
many  of  the  experimental  findings  on  these  systems.  The  bandgap  levels  for  a 
complex  of  one  to  five  chalcogen  impurities  is  shown  in  Fig.  1. 
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We  haw  da va loped  a  new  method  for  tha  efficient  calculation  of  the 
electronic  structure  in  solids.  Our  method  is  highly  physical  and  uses  a  linear 
combination  of  psaudoatomic  orbitals  in  the  first-principles  local  density 
pseudopotential  approximation.  We  have  used  this  method  to  investigate  such 
basic  properties  as  lattice  constant,  bulk  modulus,  vibrational  mode 
frequencies,  band  structures  and  charge  densitites  in  ten  different 
semiconducting  materials  of  technological  importance.  We  have  found  that  this 
tight-binding- like  technique  gives  exellent  results  when  compared  to  the  more 
rigorous  plane  wave  expansion  method.  Some  results  are  shown  in  Fig.  2. 

We  have  used  the  technique  described  above  to  study  the  formation 
energy  and  near  bandgap  levels  of  interstitial  impuritites  in  IV-IV,  III-V,  and 
II-VI  semiconductors.  We  have  discovered  an  interesting  interplay  between 
interstitial  site  preference  and  bandgap  levels  which  causes  (in  some 
circumstances)  the  impurity  to  switch  from  one  interstitial  site  to  another  as 
the  fermi-level  is  changed.  These  findings  may  have  importance  in  understanding 
diffusion  processes  in  compound  semiconductors,  and  in  determining  the 
predominant  native  defects.  The  physics  describing  this  interplay  between  the 
two  interstitial  sites  (which  we  refer  to  as  Ta  and  Tc )  is  shown  in 
Fig.  3.  Near  bandgap  interstitial  levels  for  a  number  of  defects  in  GaAs  and 
Zn&e  are  shown  in  Fig.  4 

These  calculations  have  led  us  to  work  on  predicting,  entirely  from 
first  principles,  the  equilibrium  concentrations  of  a  variety  of  defects  in 
semiconductors  including  native  defects  such  as  antisites,  vacancies,  and 
interstitials.  This  exciting  area  is  being  pursued  presently. 

The  type  of  calculations  we  are  presently  using  are  very  similar  in 
spirit  to  those  of  the  tight-binding  approximation  used  in  a  wide  varisty  of 
applications  in  semiconductors.  Our  technique  however,  is  totally  ab-initio  and 
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requires  no  fitting  to  experiment  and  has  no  adjustable  parameters.  This  makes 
it  ideal  for  the  study  of  systems  where  little  is  known  experimentally. 

Osing  an  empirical  tight-binding  method  (which  does  require 
fitting  to  experiment),  we  have  also  performed  the  first  molecular  dynamics 
simulation  of  a  compound  semiconductor  surface  using  Hellmann-Feynman  forces. 
The  calculations  use  a  parameterized  form  of  the  electron  ion  interaction,  but 
have  the  advantage  of  being  relatively  simple  to  perform.  A  unique  feature  of 
these  calculations  is  that  they  include  the  forces  generated  by  the  microscopic 
electronic  states  and  include  a  number  of  many  body  effects.  Eventually  we 
would  like  to  meld  this  type  of  molecular  dynamics  simulation  with  the  above 
ab-initio  "tight -binding"  technique.  The  relaxation  angle  $  for  the  (110) 
surface  derived  from  such  a  simulation  are  shown  in  Fig.  5.  These  results  are 
obtained  by  following  the  motion  of  the  atoms  in  time  when  a  fictitious  damping 
is  introduced  so  that  the  atoms  approach  their  eqiulibrium  configuration. 


j  We  have  also  worked  on  a  class  of  unusual  molecular  solid,  the  hydrogen 

r  halides  HF,  HC1  and  HBr.  These  materials  are  lattice  dynamically  non-linear, 


and  we  have  investigated  these  materials  theoretically  and  have  studied  their 


Fig.  i. 

Theoretical  curves  of  the  energy  levels  ( in  the  gap  between  Ev  and 
Be)  for  “tetrahedral"  nearest-neighbor  substitutional  clusters  of  from  one 
to  five  chalcogen  (S,  Se,  or  Te)  impuritites  in  Si.  The  central  atoa  in  the 
cluster  was  chosen  to  be  S  and  there  are  n  ( n=  2,  3,  4,  or  5)  impurity  atoms  in 
the  defect  cluster  around  the  central  atoa.  We  fix  the  potential  of  the  central 
atom  Vs(l)  to  that  of  sulfur,  and  vary  the  defect  potentials  Vs(x)  of 
the  other  (one  to  four)  atoms  simultaneously.  The  impurity  potential  Vs(x) 
is  negative  for  atoms  more  electronegative  than  the  host  Si  atom,  and  positive 
for  more  electroposititve  impurities.  It  is  zero  when  the  n  "defects" 
surrounding  the  central  sulfur  impurity  are  host  Si  atoms  —  in  this  case  the 
only  level  in  the  gap  is  that  due  to  the  isolated  sulfur  defect.  The  energy 
level  of  the  isolated  sulfur  impurity,  Si°,  is  driven  down  toward  the 
valence  band  when  it  is  clustered  with  other  electronegative  atoms.  Other 
states  are  also  formed  due  to  the  clustering.  The  uppermost  levels  of  sulfur 
cluster  complexes  of  2  to  5  atoms  are  shown  explicitly  at  the  defect  potential 
of  about  -5eV  appropriate  to  sulfur.  Note  the  upward  trend  toward  the 
conduction  band  of  the  highest  levels  indicating  a  continuous  transition  from 
"deep"  to  "shallow- like"  behavior  as  the  cluster  size  increases  from  2  to  5 
defects. 

Fig.  2. 

(a)  The  transverse  optic  mode  frequency  computed  entirely  from  first  principles 
using  an  sp*d*  pseudoatomic  orbitals  in  a  tight-binding-like 

hamiltonian,  compared  to  experimental  values.  The  trends  in  these  ten  materials 
are  predicted  perfectly,  and  the  agreement  with  experiment  on  an  absolute  scale 
is  excellent. 

(b)  Comparison  of  the  predicted  lattice  constants  of  semiconducting  materials 
using  sp* d*  pseudoatomic  orbitals  with  experiment. 

(e)  Same  as  (b)  but  for  bulk  modulus.  These  results  show  that  the  technique  is 


accurate  enough  to  ba  useful  in  a  variety  of  electronic  structure/total  energy 
calculations,  and  are  easy  to  implement.  The  technique  is  only  slightly  less 
accurate  than  the  more  conventional,  but  more  rigorous,  planewave  expansion 
techniques . 

Fig.  3. 

(a)  The  geometry  of  the  tetrahedral  interstitial  site,  labeled  1,  in  a 
compound  semiconductor.  The  interstitial  is  surrounded  by  four  nearest 
neighbors  labeled  1  through  4.  Two  such  sites  exist;  the  Ta  and  Tc 
interstitial  sites.  The  Ta  (Tc)  site  has  a  nearest-neighbor 

environment  of  four  anions  (cations).  One  site  ought  to  be  preferred  over  the 
other  for  a  specific  impurity. 

(b)  A  simple  schematic  three  level  diagram  for  an  interstitial  (chosen  to  be  Zn 
in  GaAs  for  discussion)  at  the  Ta  (right  panel)  and  Tc  (left  panel) 
tetrahedral  sites.  The  Zn  s-orbitals  interacts  with  the  valence  and  conduction 
bands  with  different  strengths  due  to  the  different  environments.  The  single 
dashed  line  indicates  a  weak  interaction,  while  the  double  dashed  line  a  strong 
interaction.  The  Zn  s-orbital  pushes  down  on  the  valence  level  and  up  on  the 
conduction  level.  The  circles  indicate  electrons  and  the  diagram  indicates  the 
Zn<>  center.  Note  that  the  Zn  related  "deep”  level  is  always  lower  in 

energy  at  the  Tc-site  than  at  the  TA-site,  but  that  the  change  in 
energy  of  the  valence  bands  favors  the  TA-site.  These  two  energy  terms  are 
competitive,  and  the  lowest  energy  configuration  is  expected  to  depend  on  the 
occupancy  of  the  deep  level. 

(c)  The  total  energy  for  the  Ta  and  Tc  sites  of  interstitial  Zn  in 

GaAs  as  a  function  of  the  chemical  potential  in  the  gap.  Note  that  the  bandgap 
deep  energy  levels  (indicated  by  short  vertieal  lines)  are  lower  in  energy  at 
the  Tc-site  than  the  Ta.  Also  note  that  when  the  deep  levels  are 
empty  (+2  oharge  state),  the  system  prefers  the  Ta -interstitial  site,  but 
that  the  Ta  and  Te  sites  are  nearly  degenerate  in  energy  when  the 


dMp  levels  are  occupied .  These  f natures  are  in  mgrmmmm at  with  tha  quail tativa 
physics  of  tha  diagram  in  (b). 

Fig.  4. 

(a)  Bandgap  anargy  1 avals  for  iatarstitial  impurities  ia  OaAs  at  tha  two 
tetrahedral  intarstitial  sitaa  (T*  and  Tc ) .  A  number  of  impuritias  hava  baan 
considarad  ranging  from  tha  group  I  inpurity  K,  to  tha  group  VI  inpurity  Be. 
Notice  that  tha  bandgap  is  ralativaly  "clean"  of  daap  levels  except  for  Zn,  As 
and  Sa.  Also  notica  that  tha  Tc~site  inpur i ties  hava  levels  lower  in 

energy  than  tha  Ts-site  inpur i ties.  Thasa  calculations  help  us  to 
understand  why  intarstitial  inpurities  are  not  energetically  favored  in  III-V 
materials,  but  may  be  more  energetically  favorable  in  II-VI  materials.  They 
also  show  why  tha  chemical  potential  is  quite  important  in  determining  the 
energy  of  such  defects  since  some  defects  have  a  number  of  charge  states. 

(b)  Bandgap  energy  levels  for  tha  native  intarstitial  inpurities  Zn  and  Se  in 
ZnSe  at  tha  T*  and  Tc  intarstitial  sites. 

Fig.  5. 

Relaxation  angle  4  as  a  function  of  tine  at  tha  (110)  surface  for  twelve 
seal conductors  using  a  fictitious  damping  to  determine  tha  final  equilibrium 
configuration.  Tha  atoms  ware  allowed  to  nova  according  to  Newton's  laws,  where 
tha  forces  ware  determined  by  including  changes  in  tha  microscopic  electronic 
states  via  tha  Hellnann- Feynman  theorem.  Tha  angle  Ms  a  one  measure  of  the 
relaxation  of  the  surface  atoms  from  their  bulk  atomic  posititons.  The  final 
relaxation  angle  is  found  to  be  nearly  the  same  for  all  these  materials  ranging 
from  about  25  to  90  degrees. 
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